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The recently introduced hierarchical regular networks HNS and HN4 are analyzed in detail. We 
use renormalization group arguments to show that HNS, a 3-regular planar graph, has a diameter 
growing as \/N with the system size, and random walks on HNS exhibit super-diffusion with an 
anomalous exponent — 2 — log2(0) ~ 1.S06, where (f> — (\/5 + l) /2 = 1.618 ... is the "golden 
ratio." In contrast, HN4, a non-planar 4-regular graph, has a diameter that grows slower than any 
power of A'^, yet, fast than any power of InA^. In an annealed approximation we can show that 
diffusive transport on HN4 occurs ballistically [d^ — I). Walkers on both graphs possess a first- 
return probability with a power law tail characterized by an exponent jj, = 2 — l/d^- It is shown 
explicitly that recurrence properties on HNS depend on the starting site. 

PACS numbers: 05.40.-a. 
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I. INTRODUCTION 

Networks with a sufficiently intricate structure to ex- 
hibit not-trivial properties for statistical models but suf- 
ficiently simple to reveal analytical insights are few and 
far between. Familiar examples are random graphs 
the hierarchical lattices [3] originating in the Migdal- 
Kadanoff bond- moving scheme 0, [1] , fractal lattices [1] , 
or scale- free networks 0, H • Scale- free networks and 
random graphs can elucidate mean-field properties only, 
whereas hierarchical lattices often provide excellent re- 
sults for complex statistical models in low dimensions 
but do not possess a mean- field limit. Fractal lat- 
tices provide tantalizing access to dynamical systems in 
non-integer dimensions but can not be tuned across di- 
mensions. From the perspective of statistical physics, it 
could be desirable to have a network that combines solv- 
ability for low-dimensional systems with mean-field prop- 
erties in such a way that one could interpolate between 
either extreme. 

We have recently introduced a set of networks that 
overlay a lattice backbone with regular long-range links 
[ill [T^ l , similar to small- world networks [1^] but hierar- 
chical and without randomness. These networks have a 
recursive construction but retain a fixed, regular degree. 
The hierarchical sequence of long-distance links occur in 
a pattern reminiscent of the tower-of-hanoi sequence . 
Therefore, we have dubbed them "Hanoi-Networks" and 
abbreviated them as HNS and HN4, since one is 3-regular 
and the other 4-regular. While almost identical, both 
types of networks lead to very distinct behaviors, as re- 
vealed by our studies here. Future work will focus on 
detailed studies of Ising models on these networks. In 
Ref. [ll| we have shown already that spin models on HN4 
have the desired properties mentioned above. Here, we 
analyze in great detail diffusive transport on these net- 



works, for which especially HNS proves to possess a rich 
behavior. 

The paper is structured as follows. In the next sec- 
tion, we introduce the networks and discuss various key 
design aspects. In Sec. IIII[ we explore geometric aspects 
of the networks such as their diameter. Sec. [IV] con- 
tains our analysis of random walks on both networks, 
starting with a review of the dynamic renormalization 
group in the context of a simple one-dimensional random 
walk which facilitates an efficient discussion of the equiv- 
alent but far more involved calculation for HNS. Unfor- 
tunately, the same approach does not seem to apply to 
HN4, so we conclude the section with a derivation of a 
moment-generating equation in Fourier space for walks 
on HN4 directly from the master equation and an an- 
nealed approximation. We conclude this paper in Sec.lVl 
indicating various future projects that can be developed 
from the work presented here. 



II. NETWORK DESIGN 

Both networks we are discussing in this paper consist 
of a simple geometric backbone, a one-dimensional line 
of iV = 2'^ sites, either infinite (-2'= < n < 2*^, fc ^ oo) , 
semi-infinite (O < n < 2*^, fc — > oo), or closed into a ring. 
Each site on the one-dimensional lattice backbone is con- 
nected to its nearest neighbor. (In general, the following 
procedure also works with a higher-dimensional lattice.) 

To generate the small- world hierarchy in these graphs, 
consider parameterizing any integer n (except for zero) 
uniquely in terms of two other integers (i, j), i > 0, via 



T {2j + 1) 



(1) 
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Here, i denotes the level in the hierarchy whereas j la- 
bels consecutive sites within each hierarchy. For instance, 
i = refers to all odd integers, i = 1 to all integers once 
divisible by 2 (i. e., 2, 6, 10,...), and so on. In these 
networks, aside from the backbone, each site is also con- 
nected with (one or both) of its nearest neighbors within 




the hierarchy. For example, we obtain the 3-regular net- 
work HNS (best done on a semi-infinite hne) by connect- 
ing first all nearest neighbors along the backbone, but in 
addition also 1 to 3, 5 to 7, 9 to 11, etc, for i = 0, next 2 
to 6, 10 to 14, etc, for i = 1, and 4 to 12, 20 to 28, etc, for 
i = 2, and so on, as depicted in Fig. [T] Correspondingly, 
HN4 is obtained in the same manner, but connecting to 
both neighbors in the hierarchy, i. e., 1 to 3, 3 to 5, 5 to 7, 
etc, for z = 0, 2 to 6, 6 to 10, etc, for i = 1, and so forth. 
For this network it is clearly preferable to extend the line 
to ~oo < n < oo and also connect -1 to 1 , -2 to 2, etc, 
as well as all negative integers in the above pattern, see 
Fig. [21 The site with index zero, not being covered by 
Eq. Il]), is clearly a special place, either on the bound- 
ary of the HNS or in the center of the HN4. It is easy 
to generalize these graphs, for instance, by putting the 
structure on a ring with periodic boundary conditions, 
which may require special treatment of the highest level 
in the hierarchy on finite rings. 

Random walks on these networks have fascinating 
properties due to their fractal nature and their long-range 
links as shown, for instance, in Fig. [3] 



III. NETWORK GEOMETRY 
A. Distance Measure on HNS 

For HNS, it is simple to determine geometric proper- 
ties, for instance, its diameter d, which is the longest of 
the shortest paths between any two sites, to wit, on a 
finite graph of size N = Nk = 2*-' for fc ^ oo . Clearly, 
d in this case would be the end-to-end distance between 
sites 71 = and n ^ N with the smallest number of 
hops. Using a sequence of networks for k = 2,4,6, .. ., 
the diameter-path looks like a Koch curve, see Fig. H) 
We can define the path H^ as a sequence of jumps reach- 
ing from one end of the = 2'^-long graph to the other 



FIG. 2: Depiction of HN4 on an infinite line. The center 
site is n = 0, which requires special treatment. The entire 
network becomes 4-regular with a self- loop at n = 0. Note 
that IIN4 is non-planar. 




FIG. 3: Plot of the probability Pi^t of a random walker 
to reside at a site I at time t after starting at the origin 
/ — obtained by numerical simulations on HNS. The scale 
of the vertical axis refers to In Pi_t to increase the visibil- 
ity. Note the highly non- homogeneous evolution, which is 
tied to fractal distance-dependence also apparent in Fig. [5] 
due to the hierarchy of long-range links. Those links al- 
low for walkers to appear suddenly at certain sites such as 
« = 3 X 2' = 6, 12, 24, 48, 96, 192, . . . long before ordinary dif- 
fusion along the backbone itself would provide for. Those 
onset-points for the most rapidly progressing walker should 
follow an envelope function that is linear ("ballistic") in a 
space-time plot, as walkers diffuse ordinarily with hopping- 
distance, d ~ but by Eq. (|4]), d ~ hence, I ~ t. 
Average walkers progress slower than ballistic, though, but 
still faster than diffusion. 



via 

Ho = 1, 

= l-2-l = no-2-no, 
n4 = 1-2-1-8-1-2-1 = Hz - 8- n2, 

(2) 

\ik+2 = nfc-2^+i-nfe, 

with an obvious notation of using "-" to string together a 
sequence of ever more complex moves. (While somewhat 
redundant here, introducing this notation will prove use- 
ful for HN4 below.) Hence, the length dk of each path 
life is given by 

4+2 = 2dk + 1 for Nk+2 - 4iVfc, (3) 

thus, 

d - \/iV. (4) 

This property is demonstrated also in Fig. [5l In some 
ways, this property is reminiscent of a square-lattice con- 
sisting of N lattice sites. The diameter (^diagonal) of 
this square is also ~ \/N . In this sense, presuming that 
each link corresponds to a unit distance, HNS is a fractal 
lattice of dimension 2, i. e. filling the plane. As shown in 
Fig. [51 we also find that the number of sites that can be 
reached from a given site grows quadratically with the 
number of jumps allowed. It should be noted, though, 
that we are employing the one-dimensional lattice back- 
bone as our metric to measure distances in the random 
walks in Sec. IIVI which enforces a trivial fractal dimen- 
sion of d/ = 1. We conclude that, while interesting in its 
own right, HN3 is far from any mean-field behavior for 
which we would expect typical distances to depend only 
on (some power of) InA'^. 

B. Distance Measure on HN4 

The situation is more interesting for the HN4. Using 
the notation from Eq. ([3]), we have 
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and so on. Due to degeneracies at each level (which we 
have not listed), one has to proceed to many more levels 
in the hierarchy to discern the relevant pattern. In fact. 
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FIG. 4: Sequence of shortest end-to-end paths (=diameter, 
thick black lines) for HNS of size = 2*=, k = 2,4,8. (Note 
that every second step in the hierarchical development has 
been omitted. At level k = 3, 5, . . . the shortest path is the 
same as at fc — 1, except each linear segment counts for two 
steps.) Whenever the system size A'^ increases by a factor of 4, 
the diameter d increases by a factor of ~ 2, leading to Eq. Q. 
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FIG. 5: Plot of the shortest path length between the origin 
of HNS and the Ith site on two networks of extend A'' = S 2^^ 
and = 32^^. In both sets of data, we plot the path-distance 
relative to the root of the separation between site I and the 
origin {I = 0) along the linear backbone. Then, all rescaled 
distances fluctuate around a constant mean. Those fluctua- 
tions are very fractal, their self-similarity becoming apparent 
when super-imposing the data for both sizes A" on a relative 
distance scale with l/N . 
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FIG. 6: Plot of the number of sites (or "neighborhood") Sd 
that can be reached from a given site within less than d jumps 
on HNS. Averaged over many starting sites, Sd/d? slowly con- 
verges to a constant, demonstrating that Sd grows quadrati- 
cally with d. Note that some features due to the hierarchical 
structure remain even after averaging over sites, such as the 
peaks at d = 192, 384, 768, etc. 
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FIG. 7: Plot of the shortest end-to-end path length for HN4 
networks of increasing backbone sizes N = 2^ , obtained by 
simulation. Note the piecewise-linear shape of the graph, 
which is reflected in Eq. ([8}, for example. In turn, Eq. Q 
concerns only the "bends" between each consecutive linear 
segment g. 
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any pattern evolves for an increasing number of levels 
before it gets taken over by a new one, with two patterns 
creating degeneracies at the crossover. Finally, we get 
(setting the degeneracies aside) 



Ilk = 



and so on. 

Note that the paths here do not search out the longest 
possible jump, as in Eq. ([3]). Instead, the paths reach 
quickly to some intermediate level and follow consecutive 
jumps at that level before trailing off in the end. As we 
will see repeatedly, this is the main distinguishing feature 
discriminating between HNS and HN4: Once a level is 
reached in the HN4, the entire graph can be traversed at 
that level, while in the HN3 any transport must climb 
down to lower levels (or merely jump back on that level), 
see Figs. [T]12 

Corresponding to Eq. ([6]), we obtain for the end-to-end 
shortest paths (=diameters dk) 




(fc = 2), 
(2 < fc < 5), 
(5 < fc < 9), 



(7) 



which we can generalize into a single statement introduc- 
ing a "generation" index 5 > 2, 



dl = 24:^ + (2^-1-1), 
defining h = I, where in general 



Ig-l +9, 



{lg-l<k<lg), (8) 



(9) 



demarcates the crossover point between the generations, 
see Fig. [T] Eq. easily yields 



g = 23(3 + l)-l (5>2). 



(10) 



To obtain the asymptotic behavior for d^, instead of solv- 
ing Eq. ([SI) for all fc, we note that exactly on the crossover 



points k ^ Ig (i. 



Ig-i) we have 



Defining eg = 2^d? , we get 



0-9 



(11) 



(12) 



which is easily summed up to give 

dl = (5-1)2^-1 



1. 



Remembering that k = Ig and, from Eq. ([T 
■\/2l^ ~ ^/2k ~ ■\/log2 N^, we finally get 



(13) 
that g ^ 



(N^oo) (14) 



for the diameter of the HN4. As we would expect that 
the diameter (or rather, the average of shortest paths) 
in a small-world graph should behave as d ~ log TV, it is 
instructive to rewrite Eq. (|14p as 



(log2 iV)" with a 



V2 



N 



1 



log2 log2 N 2 



(15) 



Technically, of course, a diverges with N and the diam- 
eter grows faster than any power of log2 N but less than 
any however-small power of A'^, unlike Eq. In reality, 
though, a varies only very slowly with N, ranging merely 
from a w 1.5 to 3 over fifty decades, see Fig. [5] 
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FIG. 8: Plot of the system-size dependence of the exponent 
a = a[N) defined in Eq. (|15|l . The solid curve is the exact 
value based on Eq. (|13|l . and the dashed curve is the asymp- 
totic approximation also given in Eq. (|15p . This approxima- 
tion is always an upper bound. Note that q(A'^) barely varies 
within a factor of 2 over 50 orders of magnitude. 



This is synonymous with the asymptotic form with the 
mean-square displacement 

m^t''"-, (17) 

which defines the anomalous dimension of the walk in 
terms of the exponent ■ While the result for the one- 
dimensional walk can be easily obtained with much sim- 
pler means, it serves as a pedagogical example of calcu- 
lating first-passage and first-return times using RG on 
more complex structures. Later on, reference to this pre- 
sentation will allow us to avoid excessive algebra. 

1. RG for the Id random walk 

As a pedagogical example, we present here the theory 
as it will be applied in Sec. IIV A 21 We consider a biased 
random walk on a finite one-dimcnsional line. The master 
equations for such a random walk on a lattice of length 
N = 2^+^ with reflecting boundaries arc given by 



IV. RANDOM WALKS 

In the following we will study random walks on HNS 
and IIN4. For simplicity, wc have focused in our simula- 
tions only on two elementary observables, the mean dis- 
placement with time, {\l\) ^ t^/d^ ^ and the first-return 
time distribution Q{At) ^ Ai~''. All walks are con- 
trolled by one parameter, p, which determines the proba- 
bility of a walker to step off the lattice into the direction 
of a long-range jump. In particular, on HNS a walker 
steps off the lattice with probability p and jumps either 
to the left or right neighbor with probability (1 — p) /2, 
whereas on HN4 long-range jumps to either the left or 
right occur with probabihty p/2 instead. In both cases, 
we should return to a simple one-dimensional nearest- 
neighbor walk for p —^ 0, where d^, = 2 and fi = S/2, al- 
though we would expect that limit to be singular. Here, 
we only consider uniform values of p, independent of the 
sites (or level of hierarchy). With a bit more algebra, the 
following considerations could be extended to values that 
for each site depend of the level of hierarchy, p = p{i), 
for instance. 



A. Renormalization Group for Random Walks 

Here, we analyze the spatio-temporal rescaling of sim- 
ple random walks with nearest-neighbor jumps along the 
available links using the renormalization group. First, we 
review the method using as a simple example the well- 
known one-dimensional walk. Walks on such a graph 
obey simple diffusion, ~ 2, which implies that each 
rescaling of space entails a rescaling of time according to 

N ^ N' = 2N T^T' ^ 2'^^T = 4T. (16) 



Po.t+i = pPi,t, 
Pi,t+i = Pa.t+pP2,t, 

Pu+i = {l-p)Pi-u+pPi+i,t i2<l<N-2), 

PN-l,t+l = {l-p)PN-2,t+ PN,t, 

PN,t+l = {l-p)PN-l,t, (18) 

where Pi^t denotes the probability of a walker to be at site 
I at time p is the probability expressing the biasing for 
left or right hops. Since we want to start the walks at 
time t = at the origin / = 0, these equations have the 
initial condition 

Pi.Q = 5i,a. (19) 

To facilitate renormalization this non-equilibrium pro- 
cess, we introduce a generating function 

oo 

Pi{z) = Jl^i.t^* (20) 

for all < ^ < A''. Incorporating the initial condition in 
Eq. (Ull), Eqs. HH) transform into 

Pi = cP^ + aP2, 

Pi = hPi^i+aPi+i {2<l<N-2), 

Pn-1 bPN-2 + dPN, 

Pn = bpN-i, (21) 

where we have inserted generalized "hopping rates" in 
preparation for the RG. Initially, at the k = 0th RG 
step, it is 

a(°) = pz, 

= {l-p)z, (22) 
c(°) = z, 

= z. 
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which provides a sufRcicnt number of renormahzable pa- 
rameters that are potentially required to consider special 
sites at both boundaries. 

A single step of applying the RG consists of solving 
Eqs. (|2ip for Pi with odd values of / (which is trivial 
here, as they are already expressed explicitly in terms of 
even ones) and eliminating them from the equations for 
the even I. After that elimination, we can rewrite the 
equations for even I as 



Po 

P2I 

Pn~2 

Pn 



■P2 



1 



1 — ac 1 — ac 



be 



P 



l-2ab 1 - 2ab 



P21- 



■P4 



1 - 2ab 



1 - 2ab 



P21 



+2 



TV 

2<l<- 



&2 ad ~ 

PN-i + -, — :r-rPN, 



l-2ab'" " ' l-2ab 
62 



1-bd 



P 



N-2- 



Comparing these equations with Eqs. (PT|) allows to ex- 
tract the RG recursion equations. [Note that superscripts 
referring the A:th RG step have been suppressed thus far 
in Eqs. ^ and 

Before we analyze the first return time at the bound- 
ary specifically, we can use the equation for bulk sites 
I in (|23p to extract already the diffusion exponent dw 
A comparison of the respective coefficients in Eqs. (PT|) 
and ((^ yields 



^(fc+i) 



^(fe+i) 



l-2d^)b^^)' 

V 2 



(24) 



y 



(fc) 



y 



(k) 

Vi e 



(25) 



Inserting the Ansatz in Eq. ([^5]) into the recursions in 
Eqs. ([24]), we obtain near the fixed point with a* = b* = 



I. 
2 ' 



(fe+i) 

dk+l) 



(fc) 



ib\ 



(fc) 



(26) 



with the result that 



6^'■)(x4^ 



This implies that as space rescales by a factor of 2 (i. e., 
eliminating all odd-index sites), time rescales by a factor 
of 4, as indicated in Eq. (|16p for an unbiased random 
walk, leading io d^ =2. The same analysis for either of 
the biased fixed points yields that, for example, a' 
beyond any power of e and, with the Ansatz 6'^'^^ ^ 1 







b'^^^e in Eqs. ([M]) . b'^^'' oa 2^ . Following the interpretation 



in Eq. (jl6[) . this would imply d^; = 1 and we find the 
famihar result that with the sHghtest bias, i. e., p < ^ 
or p > i, the motion at large length and time scales is 
\dominated by the constant-velocity drift upon reaching 
2 Jthe bulk. 

In this scenario of a bias, average first-return times are 
clearly system-size independent constants: A walker with 
a bias towards the origin {p > ^) will drift back recur- 
^2j^iitly after only small excursions; a walker with a bias 
away from the origin (p < i) returns at most a finite 
number of times in short order until the drift eventually 
carries it away without further recurrence. In the follow- 
ing, we therefore focus exclusively on the unbiased case 
p = ^. Then, we can equate a = & at every step, to get 
from Eq. 



These recursions converge for fc 00 towards fixed 
points (a*, 6*) that characterize the dynamics in the 
infinite-time limit (which corresponds to the limit of 
z — > 1~). The trivial fixed point a* = 6* = is unphys- 
ical, as it can not be reached from the initial conditions 
in (22) for any choice of p (and z = 1). The physical 
fixed points are (a*, b*) ~ (1,0), which is reached for any 
bias p > i, or (a*, b*) = (0, 1), reached for p < 5; finally, 
(a*, 6*) = (^,^) can only be reached by entirely unbi- 
ased walks for p = i. To explore the behavior for large 
but finite times, we expand the RG recursions in ((24)) to 
first order in e = 1 — z by writing for y G {a,b}: 



l-2{aik)y 



(27) 



To derive the return-time behavior, we have to examine 
Eqs. ([23]) more closely. Comparing coefficients also in the 
boundary terms leads to 



rf(fe+i) ^ 



a 



l-2(aW)^ 
l-2(aW)' 



(28) 



For large fc, both S'^^ and d*^*^^ are entrained with a^''\ 
and we obtain a consistent and closed set of relations 
for all coefficients in Eqs. (^J) and (221) by identifying 
c = d ^ 2a. Further renormalizing 



p(fc+i) 



l-2(a^^^)' 



p(k} 
^21 



(29) 



ensures invariance of the constant term at the lower 
boundary that originated from the unit initial condition 
in Eq. (fTOll . 

After fc — K RG steps, the system has reduced to 



JK)p(K) 



P. 



2a^'^>P^'^' +2a'-'^>P^ 



(30) 



7 



which yields 



Pr 



l-2(aW)' 
1 - 4faWV 



(31) 



From the ratio of {t) and (i^)^, we obtain then 



/X = hmjv-,oo 2 



(38) 



Using Eq. (P^j) in turn obtains 



p(0) _ 



n 



fe=0 



1 



l-4(a(^))^ H [l 



(32) 



It is a well-known fact that the generating functions for 
being at the origin, and for the first-return proba- 

bility to the same site, Qoj satisfy the following simple 
relation [Til: 



= 1 - 



1 



(0)' 



(33) 



Note that a recurrent walk (with Qo = 1) requires that 

^1^). In our one- 
* = i , on behalf of 



Pq"' diverges at long times (i. e. z 
dimensional walk here, it is a(^) ^ a* 

which the denominator of P^^ in Eq 
making the walk recurrent. In more detail, it is 



has a zero, 



Qo = 1 ^ — '-2 TT 

l-2(aW)^ H 



1 - 2 ( aC^) 



(34) 



We intend to extract the exponent fi for the first-return 
probability distribution, which on a finite but large sys- 
tem of size N oo behaves as 



(35) 



where tjv is a cut-off time scale that diverges in some 
form with N. As the physics of the return probabilities 
at any finite time should not change if the system size 
becomes infinite independently, it must be /i > 1 for 
Qajt) to remain normalizable. Based on that observation 
jl4| . we need to calculate the first two moments of Qo{t), 
corresponding to an expansion of Qo{z) to 2nd order in 
e = 1 — z, to extract /i. Since /i > 1, the normalization 
integral 



Luckily, due to the leading zero in 1 — 4 (a^-'^^) , any 
other factor in Eq. ([M]) only needs to be expanded to 
first order in e. Extending the Ansatz in Eq. for a^*^' 
to 2nd order, we obtain here: 



i - i X 4^'^e + 5 X 



K2 



1-4 



A' 2 



1 - 2 ( 



K-l r 



n 

fc=0 



1 - 2 a 



X 4^6 + 116 X 16^6 



- + 4 X 4^e -f . . . , 



n 

k=0 



4 X i^e 



K-l 



1 + 8e ^ 4* 



fc=0 



2^e 



Inserting these expressions into Eq. (|34p . we get 



1 - 16 X 2^e 



952 



X 8^^e2 ^ 



(39) 



The moments of Qoit) are obtained via derivatives of 
Qo{z), i. e. 



dz 



(40) 



Applied to Eq. ([39|) . we calculate for the dominant 
asymptotics of the moments: 



2^ - N, 



which from Eq. (|38p leads to the familiar first-return ex- 
ponent of a one-dimensional walk. 



(41) 



dtr 



0(1) 



(36) 



is dominated by its behavior for small t, which is irrele- 
vant in detail except for the fact that it makes the integral 
become a constant independent of t^t. If we further as- 
sume that < 2, then for all m > 1, the integrals for 
those mth moments do diverge with r^r. 



IN 



1 

AT 



(37) 



2. RG for the random walk on HN3 

We follow the discussion for one-dimensional walkers 
above to model diffusion on HNS. The lesson of the pre- 
vious section is that we only need to consider the bulk 
equations to extract the diffusion exponent du,, and sup- 
plement with the calculation on the final system at the 
end of the k = Kth step of the RG to obtain the first- 
return exponent for a system of size N = 2^+^. We can 
section the line into segments centered around sites I with 
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i = 1 in Eq. i. e. n = 2(2j + 1). Such a site / = n is 
surrounded by two sites of odd index, which are mutually 
hnked. Furthermore, n is hnked by a long-distance jump 
to a site also of type z = lat/ = n±4in the neighboring 
segment, where the direction does not matter here. The 
sites I = n±2, which are shared at the boundary between 
adjacent segments also have even index, but their value 
of i > 2 is undetermined and irrelevant for the immediate 
RG step, as they have a long-distance jump to some sites 
I = m± at least eight sites away. 

In each segment in the bulk, the master-equation reads 



Pn+l.t + 1 
Pn,t+1 
Pn-l,t + l 
Pn-2,t + l 



1 




P 




2 




1 




P 




2 




1 




P 




2 




1 




P 




2 




1 




P 


2 



[Pn+3,t + Pn+l,t] + pPm+,t, 
[Pn+2,t + Pn,t] + P Pn-l,t, 
[Pn+l,t + Pn-l,t] +pPn±4,t, 
[Pn,t + Pn-2,t] + P Pn+l,t, 
[Pn-Ut + Pn-3j] +pPm^,t- 



(42) 



Here, p is the (uniform) probability for a walker to take a 
long-range jump, whereas (1 — p)/2 is the probability to 
jump either left or right towards a nearest-neighbor site 
along the backbone. Without restriction of generality, let 
us assume that we happened to let all walks start from 
a site within this segment. Note that, unlike for the Id- 
walk above, there are three distinct types of sites even 
in the bulk of this problem: Pn±2,ti Pn±i.ti and Pn,t- 
In principle, one could expect three distinct return-time 
behaviors as a result. We will demonstrate below that, 
indeed, different sites possess differences in their return- 
time behavior, depending on their level in the hierarchy, 
but all scale with the same exponent /z. Here, we choose 
to have the walk start on site n — 2 within this segment: 



Pi.o 



(43) 



Using the generating function in Eq. (PO)) on Eqs. (1421 - 
yields: 



Pn+2 


~ a 


Pn+3 


+ Pn 


fl 


+ c 


Pn^ 


h4 + Pn 


+ P2 Pni+ , 


Pn+1 


= b 


Pn+2 


+ Pn 


4 


Pi Pn-1, 






Pn 


— a 


Pn+1 


+ Pn 


-1 


+ c 


Pn- 


V2 + Pn- 


-2 


+ P2Pn±i, 


Pn-1 


= b 


Pn + 


Pn-2 


4 


Pi Pn+1, 




(44) 


Pn-2 


= a 


Pn-1 


+ Pn 


-3 


+ c 


Pn 


+ Pn-4 


+ P2 An_ +1, 



where we have again absorbed the parameters p and z 
into general "hoping rates" , which are initially 



c(0) 

(0) 

Pi 



0, 

Pf 



(45) 



P2 

b 



Pi 



n-2 



n-l 



9 



9 



n+l 



9 



n+2 



c' c' 



9 



n-2 



9 9 



n+2 



FIG. 9: Depiction of the (exact) RG step for random walks 
on HNS. Hopping rates from one site to another along a 
link are labeled at the originating site. The RG step consists 
of tracing out odd-labeled variables Pn+\ in the top graph 
and expressing the renormalized rates (a', fe', c',p'i,p2) on the 
bottom in terms of the ones {a,b,c,pi,p2) from the top. The 
node Pn, bridged by a (dotted) link between P„-i and Pn+i, 
is special as it must have n = 2(2j -I- 1) and is to be decimated 
at the following RG step, justifying the designation of p'l. 
Note that the original graph in Fig. [T]does not have the green, 
dashed links with hopping rates (c, c'), which emerge during 
the RG recursion. 



Note that we have added new terms with a parameter 
c, which is zero initially. As Fig. [5] depicts, such links 
arc not present in the original network HN3, but have to 
be taken into account during the RG process. In all, a 
surprising number of parameters, five in all, is required 
even for the most symmetric set of initial conditions to 
obtain a closed set of RG recursion equations. Unlike 
for the one-dimensional walk in Eq. (j2ip . the parameters 
a and b here do not express a directional bias or drift 
along the backbone. Instead, they are necessary merely 
to distinguish between hops out of (currently) odd and 
even sites, respectively, as shown in Fig. [5) 

The RG update step consist of eliminating from these 
five equations those two that refer to an odd index, n± 1. 
As a first step, adding the two relations referring to the 
indices n ± 1 , we obtain 



Pn+1 + Pn-1 — 



1 -pi 



Pn+2 + 2P„ + Pn- 



zp. 



which allows us to eliminate any reference to odd sites 
n ± 1 from the middle relation in (|44p . Furthermore, 
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solving for P„±i explicitly, 

b ~ bpi ~ b 

Pn±l = 2 Pn±2 + "j J ^"T2 + ' 



I -Pi 



I -Pi 



i-Pi 



and inserting into the relations for Pn±2 in Eq. 
suits in 



re- 



n+2 



[a6 + c(l-pi)](l+pi) 



l-pj- 2ab 
abpi 
l-pl-2ab 

P2{1-pI) 



Pn+A + Pn 



Pn+6 + Pn-2 



.2 o„;,^™+' 



Pn = 



Pn-2 — 



I -pi- 2ab 
ab + c{l — pi) 



1 — pi — 2ab 

^ P2(l-Pl) 

1 — pi — 2ab 
[ab + c{l-pi)]{l+pi) 
1-pI- 2ab 
abpi 



Pn+2 + Pn-2 
Pi ±4, 

Pn + Pn-4 



(46) 



Pi 

P2{i-pi) 

l-pl-2ab 



P 



'n+2 + Pn-& 



P 
± n 



i-pi 



I-pI- 2ab' 



Given that n is only once divisible by 2, either n — 2 or 
n + 2 must be divisible by 2 at most twice, and we assume 
without restriction of generality that n + 2 satisfies this 
description, whereas n—2 is divisible by a higher power of 
2. (This also selects the upper sign in the index "n±4"). 
As in Eq. ([29|) . we now add the proper superscript for 
the fcth RG step and obtain the renormalized generating 
functions at step A; + 1 as 



(,+1) _i-^^^y-2os>^ 



p(fe) 



(47) 



The proportionality factor in this relation arises from the 
necessity to preserve the unity of the initial condition, as 
in Eq. ([29]) above. Note that we would have obtained the 
identical factor, if we had chosen to start the walker from 
the central site n instead of n + 2. Starting at a site like 
ri ± 1, in turn, would have provided a different factor. 



(fc) ' 



1 



(48) 



-Pi 



potentially leading to a distinct return-time behavior. 
Yet, asymptotically for large times and distances both 
alternatives prove identical to sufficiently high order as 
to not affect the scaling discussed below. 

Proceeding with the RG step, the role of the central 
site in the new segments is then specified via 



n+2 



n such that P, 



(k) 



P. 



(fc+i) 



1 n 



(fc) 



P 



(fc+1) p{k) 



n+2 



(fc+1) 



and correspondingly for the functions on the left- 



P 

hand side of Eqs. which now reads 



p(k+l) 



p(fe+l) I p(fc+l) 
ri+1 ' n-1 



p(fe+l) , p(fc+l) 
^n+2 + ^n-2 



P2 



(fc+1) p(fc+l) 



P 



(fc+1) 



P 



Tl-1 
(fc+1) 



n-2 



p(fe+l 



P. 



n-2 



Pi 



n±4 I 
(fe+1) p(fe+l) 



n+1 J 



p(fc+l) , p(fc-+l) 
^n-1 n-3 



(49) 



P( 



fc+1) 



P 



(fc'+i) 



+p(fc+i) /,(fc.+i) + 1^ 

These equations have exactly the desired form of the cor- 
responding unrenormalized ones in Eqs. (|44p . necessitat- 
ing renomalization recursions for the parameters of the 
form 



,(fc+i) 



5(fc-+i 



.(fc'+i) 



(fc+1) 
Pi 



a(fc)6(fc)+c(fc) (l-pf^ 



1+K 



(fc) 



j(fc)&(fc) +cW 



-Pi 



(fc) 



\ - p^^^ - 2aim^) 
a(fc)Mfc)pW 



(fc) 
Pi 



2ai^')b(^') 



(50) 



pf' {I-Pf') 
I - p^^^ - 2aim^) 

1 



(fc) 

P2 



(fc) 
Pi 



1 - [pf^^ -2aW6(fc) 



The analysis of the fixed points for fc ^ oo of Eqs. (|5D)) is 
surprisingly subtle. Of course, we obtain a rather simple 
fixed point for the choice of p = 0, which eliminates all 
long-range jumps. Then, c'*^', p''i \ and p^^p vanish for 
fc = in the initial conditions in (j45p and remain zero for 
all fc > 0, according to Eqs. ((50l) . As a consequence, the 
distinction between a'^^'^ and 6'^'^^ disappears and both 
recursions reduce exactly to that for the unbiased one- 
dimensional walk in Eq. ((27)l with a* = &* = i, leading 
to ordinary diffusion with (1^ — 2 and fj. — | , as discussed 
in Sec. IIV A II Clearly, this fixed point is unstable with 
respect to variations in p. 

For any probability p > inserted in Eqs. (j45|) . the 
recursions in (|50p evolve towards an apparent fixed point 
at a* = 6* = c* = and pi = P2 = 1- But these 
recursions are singular at such a fixed point, requiring a 
more detailed investigation. If we choose an arbitrarily 
small (5 > and set p = 1 — S, we find that for all fc > it 



i(fc) ^ 6(fc) ^ c^i^) = 0{S) and pi' 



(fc-) 



(fc) 
P2 



1-0{S). 
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Hence, setting (5 = in the end indeed validates this fixed 
point. Yet, the physics of this fixed point, corresponding 
to choosing p = 1, is trivial and does not reflect the 
numerical observations: Strictly for p = 1 there is no 
transport at all along the backbone, and any walker is 
confined forever to jump back-and-forth on the first long- 
range link it accesses, which would imply dui = oo in 
Eq. pT|) . For any choice of < p < 1, no matter how 
close to unity p gets, at long-enough times the walker 
"escapes" along the backbone sufficiently often to explore 
ever-longer jumps both, to prevent confinement and to 
exceed ordinary diffusion. Thus, we must conclude that 
even this confinement fixed point is unstable and there 
has to be a third fixed point, at least. 

To find this fixed point, we have to move beyond look- 
ing at the stationary behavior (fc = oo) of Eqs. ([50)) . Al- 
though Eqs. ([50)) represent a five-dimensional parameter 
space, there do not appear to be any further stationary 
points reachable from the initial conditions in Eqs. ([15)) 
aside from those two already discussed. As all flow ap- 
pears to converge towards the singular confinement fixed 
point, we make an Ansatz inserted into Eqs. ()50p for 
k ^ 1 that explores asymptotically the boundary layer 
[l5| in its vicinity: 



y 



(fc) 



{y e {a, b,c,l~ pi, 1 - P2}) , (51) 



with the assumption that a > 1. Expanding Eqs. ([50)) to 
leading order in a~'', we find an over-determined system 
of equations 



It has a unique solution that satisfies the condition on 
the boundary layer, a > 1, namely 



1.2361 



(56) 



Thus, we found another fixed point, lurking in the 
boundary layer of the confinement fixed point but with 
very distinct physical properties from it. Each reduc- 
tion of the system size by a factor of 2 is accompanied 
by a rescaling of the hopping parameters by a factor of 
a~^, bringing them closer to confinement, yet, leaving 
just enough room to escape and find still longer jumps. 
In fact, at this point of the analysis it is not even obvious 
whether these ever less frequent escapes from total con- 
finement ultimately would result in sub-diffusive, normal, 
or super-diffusive behavior. 

The leading-order Ansatz in Eq. (|5ip merely provided 
the existence of a third fixed point at infinite times. 
Extending the analysis to include finite-time corrections 
(i. e., e = 1 — z ^ 1), we include first-order corrections. 



y 



(fc) 



Ay a 



-k 



{l + eS,/3^- + ...}, 



(57) 



expand the recursions in ()50p in e, and then also use the 
fact that is small. Using the leading-order constants 
Ay in Eqs. ([55)) . also the next-leading constants By are 
determined self-consistently. We extract, again uniquely. 



2.^1 



(58) 



a 



-Ar 



Al 



Ar. 



—M-pt 

a 



-Al^r 



(52) 



Ai-pi 



Exercising the freedom to choose Aa = 1, vfe find 

Al, = Aa = 1, 

" 2^' 
Ai-pi = 2, 
Ai-p^ = 0^, 



(53) 



where 



75-1-1 



1.6180. 



(54) 



is the "golden ratio" [l6| , and obtain the eigenvalue equa- 
tion 



a^{a + i) = 8. 



(55) 



and find, choosing = 1, 



Bh 



B„ = 1, 



6 



Bi-pi = -7, 
B1-P2 = 

Accordingly, time rescales as 
T ^ T' 



5' 
8 

502- 



and we obtain from Eq. ()17|) with T • 
sion exponent for HN3: 



(59) 



(60) 



U^'- for the diffu- 



= 2-log2</> = 1.30576. 



(61) 



Our simulations in Fig.[TO]are in excellent agreement with 
this result for d^,. 

To extract the scaling behavior of the first-return dis- 
tribution Qo{t) defined in Eq. ([35]), we proceed similar 
to the discussion in Sec. IIV A II Here, we will find 
that different sites behave differently with respect to their 
return-time behavior, which is not surprising as the hier- 
archy of long-range jumps restricts translational invari- 
ance along the backbone. As in the discussion of first 
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FIG. 10; Plot of the results from simulations of the mean- 
square displacement of random walks on HNS displayed in 
Fig. [T] More than lO'^ walks were evolved up to fmax = lO'' 
steps to measure {r^)t- The data is extrapolated according 
to Eq. (|17|l . such that the intercept on the vertical axis de- 
termines dw asymptotically. The exact result from Eq. (|6ip 
is indicated by the arrow. 



returns in Sec. IIV A 11 we anticipate that we need an 
expansion of the parameters to order e^, i. e., we extend 
the Ansatz in Eq. (|57p even further. We find that it is 
sufficient to use 



y 



(fc) 



Aya-^ {1 + eByp'' + + . . .} 

Ay |a~'^' + eBy2^ + e^Cy {^af + . . .} (62) 



with 



Ca 

C'l-pi 

Cl-P2 



121 
25" 16' 



121 0^ 

"25" 16' 
121 

~"25" 16' 
121 1 

^"25" 16^' 



(63) 



where we have also dropped at each order in e terms 
that grow less than the leading exponential in k. For the 
coefficients Cy there is no freedom to choose, as they are 
a result of quadratic terms of the previous order, eByP''; 
other terms with that freedom are subdominant. 

As for Eq. pO|) . after k = K — 1 RG steps, the system 
has reduced to an elementary graph consisting of a single 
segment like that one shown in Fig. [5] (left) with n = 2 
but with p2 = 0. Now. all even sites {I ~ 0,2,4) are 
no longer connected to a long-range jump. (Jumps of 
rate c do not count as long-range, since c — > for k — > 
00, whereas pi ~ P2 1-) We therefore consider two 
(extreme) possibilities: (1) Returns to a starting point 
at the boundary {I = 0,N) or central site {I = A^/2) 
on the network, and (2) returns to a site {I ~ N/A or 
I = iN/A) with the longest-possible long-range jump on 
the original network. 



First, we consider case (1) of starting at a boundary 
site, say, 1 = 0. It is easy to show that the central site 
I = N/2 behaves identical to those on the boundary, 
and furthermore, even making the boundary sites more 
accessible be using periodic boundary conditions does not 
change the conclusion. We solve the system of equations 



-Pi 
P2 

Pi 



aPi+cP2^ 
2bPo + bp2 



1, 

-Pi P3, 



2c 



a [Pi + P3 

bP2+2bPi+piPi 

aP3 + cP2, 



P0 + P4 



(64) 



where we have suppressed the superscript (^~^) on the 
generating function and the parameters alike. [As we 
have learned from the Id- walk in Eq. the hopping 

parameters originating from the boundary sites / = and 
I = 4 have to be doubled due to the reflecting bound- 
aries.] The solution for Pq of this system of equations 
yields 



Po 



[1 - 2^) (1 - pI) - 2ab [1 (1 -h 2c) (1 + pi)] 



(65) 
2a2fe2 



[(1 - pi) (1 - 2c) - Aab] (1 + 2c) (1 + pi - 2ab) 

Notice that both, numerator and denominator, decay 
asymptotically like ^ 1 — p\ near the stable fixed point. 
Dividing out that behavior, both behave as 1 -t-y*-^-* with 



y 



{K) 



as in Eq. where y here stands for a mix of 

coefficient that results from multiplying out the terms in 
numerator and denominator, respectively. Since we do 
not expect any spurious cancellations, the ratio of these 
two expressions also results in the form 1 -\-y'^^\ Hence, 
making superscripts reappear, we find from Eq. (p^ 



1 + Se2^+Ce'"(4a)^ -t- 



K 



(66) 



where we have also dropped the a~^-term in order e° and 
marked constants that are unimportant for the scaling 
with K by calligraphy script. Using Eq. (|47l) inserted 
into the relation for Qq in Eq. (j33p . we find 



= 1 - 



1 



P 



(0) 



K-2 



1 - 



P 



(K-i) n 



fc=0 



1 - 



2a(fc)6(fc) 



1 - bi 



(67) 



Ignoring at most a finite number of factors in the product 
(hence, missing an overall constant A), we can expand 
the remaining factors in the product using the asymptotic 
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expansion in Eq. ([62)) for 1 ^ fc < if — 2 oo: 



K~2 

n 



fe=0 



1- Pi 



X'e2'= +fe^ (4a) 



All [l-a-^ 
fc»i 



A + Te2'' + ge' (4a)" + 



(68) 



Note that < ^ < 1, as each factor in the product 
must be between i and 1 for any choice the probabihty 
p. Inserting Eqs. ([66]) and ([68]) into Eq. ([67|l . we obtain 



go ~ (l-A)+ne2^+Ie^{Aa) 



(69) 



This result imphes first that returns to sites hke / — 
0, iV/2, N are not recurrent because walks may escape 
forever with a finite probability A. Finally, an applica- 
tion of Eq. (|40|) produces 



(t'), - (4a) 
hence, with Eq. ((38)) and a = 2/0: 

1 



2 - log2 (0) 



1.23416. 



(70) 



A relation of this form, = 2 — l/d^, is commonly found 
for Levy flights [l3l , and the result is again borne out by 
our simulations, see Fig. [TTJ 

Case (2), referring to walkers starting near the longest 
jump in the network, is represented by the system 



Po = aPi+cP2, 

Pi = 26P0 + 6P2+P1P3 + I, 



P2 = a 



2c 



P0 + P4 



(71) 



P1+P3 

P3 = bP2 + 2bPi+piPi 

P4 = aP3 + cP2, 



with the constant term from the initial conditions now 
near an odd site, say, 1 = 1. Solving for the launch site 
of the walk, we find 



Pi 



1 - 2c - a6 (3 - 2c) 



[(1 - pi) (1 - 2c) - Aab] {1+pi- 2ab) 



(72) 



Note that in this case only the denominator of Pi = 
decays with 1 — pi on approaching the stable 



Pi 



{K-l) 




FIG. 11: Plot of the probability Q{At) of first returns to the 
origin after At update steps on a system of unlimited size. 
Data was collected for three different walks on HNS with p = 
0.1 (circles), p = 0.3 (squares), and p = 0.8 (diamonds). The 
data with the smallest and largest p exhibit strong transient 
effects. The exact result in Eq. (|70|l . fi — 1.234 . . ., is indicated 
by the dashed line. 



this process recurrent on behalf of Eq. (|55)) . In contrast 
to Eq. (|66)) . we find here 



1 



PI 



(K-l) 



(73) 



(I-Pi) 1 + Je2^ + /Ce"(4a)" + 



Ai-p^a 



-K+l 



l + Je2^ + ICe^ (4a)^ + 
1 + eBi^p, i2af~^ + e^Ci^p, {2af^^-^ -I 
(2a) a-^ \l + Ce {2af +Me^ (2a)'^' + 



where the last step is justified by the observation that 
the leading terms order-by-order in e in the first square 
bracket dominate over those in the second. 

We have to apply Eq. ((48)) for this case to obtain 



K-2 



(fc) 



=0 l-pf^ -2aW&W 



and from Eq. ((3S 



Qn = 1-- 



(0) 



K-2 



■'I 



2a('=)6« 



1-K 



(k) 



(74) 



fixed point which will ensure that P?,^ diverges, making all form of the expansion when multiplying l/P\^ in 



The product in Eq. ((7^ behaves identically to that dis- 
cussed above in Eq. ((55)) and hence will not alter the over- 

{K-l) ■ 
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Eq. dZSl). Then we get 

Qn - 1- {2a)a-^ \l+J\fe {2af + Oe^ {2af'^ + . 

~ l~{2a)a-'' +Ve2"' + Qe^Aaf + .... (75) 

This result is almost identical to that for returns to the 
origin above in Eq. ((69|) . except that sites at the longest 
jump in the system are recurrent, i. e., Qjv/4 = 1 for 
large systems (K — > 00) and times (e = 1 — z 0). 
(Presumably recurrence will gradually degrade from the 
strictly recurrent sites at the highest level in the hierarchy 
to those at the lowest.) Yet, the scaling of return times 
in form of the exponent fj, is described by Eq. ([70]) for all 
sites. 



B. Generating Function for Random Walks on 
HN4 

Next, we consider a random walk on HN4. The 
"master-equation" [l3| for the probability of the walker 
to be at site n, as defined in Eq. ([1]), at time t is given 

by 

+ 1 [Pn-2^+i,t-l + P„+2'+i.t-l] , (76) 

where p is the probability to make a long-range jump. 
(Throughout, we considered p uniform, independent of n 
or t). To make the connection between n and i explicit, 
we rewrite Eq. (j76p as 



Po, 



l-p 



(2j + l)-l,t-l 



P 



2i(2i+l) + l,t-lJ 



2'(2i-l),t-l 



P 



2'{2j+3),t-l\ J 



l-p 



Pi 



t-l\ 



pPo.t 



(77) 



Note that the case n = is not covered by Eq. ([T]) and, 
hence, must be treated separately. Here, we choose the 
site n = to be the only one connected to itself such that 
HN4 is 4-regular throughout, as depicted in Fig. [21 

It is straightforward to apply the generating function 
in Eq. ((20)) again, assuming, for simplicity, the initial 
condition 



(78) 



We obtain 



^2^2^ + 1) 



Pn-1 



l-p 



1 P 

l-p 



-P2'(2j + 1)-1 



P2'{2j+1) + 1 



P 



2'(2j-l) 



P 



P-1 + Pi 



2'(2j+3) 



pzPo. 



(79) 



While the overall structure of this problem is even more 
symmetric than for HN3 in Sec. IIV A 2) a RG treatment 
does not seem possible in this case. Tracing out all odd 
sites would immediately interconnect all other remain- 
ing sites. (The resulting infinite set of coupled equations 
may have certain symmetry properties that would lend 
themselves for a recursive treatment. We have not yet 
explored such a possibility.) 

In contrast to an ordinary lattice, say, it is also not 
straightforward to solve this equation by a Fourier trans- 
form such as 



F(z,0) = J2 ^"(a:)e"^^ 



(80) 



defining / = \^—T- Considering the l-p terms, originat- 
ing from nearest-neighbor jumps, and the p terms, orig- 
inating from long-range jumps, in Eqs. (j79p separately 
provides for regularly-space patterns, level-by-level in the 
hierarchy. But the mixing of nearest-neighbor and long- 
range jumps destroys this regularity. Hence, we resort to 
transforming the Eqs. (j79p in each level i with a partial 
transform, 

00 

n,(z,0)- ^ P2 '{2]+i){z) exp {2* {2j + 1) (j)l} . (81) 

j=-oc 

Inserting Eq. (|8T|) and application of the general theorem 



00 00 



(2i+i)±i 

2—1 j — — 00 



E •^2j + l 

j=-oo 



/±i, (82) 



which results because 2*(2j -I- 1) for i > 1 exactly runs 
over all even numbers 7^ and over all odd numbers for 
i = Q, yields for Eqs. ((79)) : 

00 

^Hi [l-p0COs(2*+V)] 

= (1 - p)z COS (j)Ila + 1 - (1 -pz)Po, 

Ho [1 -pzcos(20)] (83) 

00 

= (1 — p)z cos (j> Hi + (1 — p)z cos (/)Po- 

i=l 

We can combine both relations to get 

1 = [1 — zp— {1 — p)z cos (j)]F 



~zp 



£[l-cos(2'+V)]n., (84) 



i=0 



using Eqs. ((50)) and ((5T)) to eliminate Pq via 



00 00 



En^ = E E ^'2=(2j + l) cxp{2^(2j + l)0/} 
i=0 i—0 j — — oo 



71— — 00 

F-Po. 



(85) 
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At this point, there does not seem to be any fur- 
ther progress possible on Eq. fM]) . due to the term 
Ilicos (2*+^0) , which resembles a Weierstrass func- 
tion [l3l- At best, on could try to extract information 
about the moments of the walk, 

oo 

via the moment-generating function 

oo 

M,(z) = Y.(^')t''^ 

t=Q 

= [-/9^]'=F(z,0)U^o. (86) 

Note that the 2nd moment Al2{z) already would provide 
the exponent dw on behalf of the definition in Eq. (fT7|) . 
[All odd moments vanish, of course, as Eq. ([84]) is even 
in (j).] The 0th moment, setting = in Eq. simply 
results in 

Moiz) = 

1 — z 

which just demonstrates that everything is properly nor- 
malized, J\ft = = J2n=-oo Pn,t = 1, at all times t. 
But already the 2nd moment would lead to terms con- 
taining ni4% which we can not account for, even at 
(/) = and in the limit z ^ 1 ^ . 

Instead, we note that the long-time behavior is dom- 
inated by the long-range jumps, as discussed for HNS 
in Sec. IIV A 21 To simplify matters, we set p = 1/2 
here, although any finite probability would lead to the 
same conclusions. We make an "annealed" approxima- 
tion, i. e., we assume that we happen to be at some site 
n in Eq. U]) with probability 1/2', corresponding to the 
relative frequency of such a site, yet independent of time 
or history. This ignores the fact that in the network ge- 
ometry a long jump of length 2' can be followed only by 
another jump of that length or a jump of unit length, 
and that many intervening steps are necessary to make a 
jump of length 2'+^, for instance. Here, at each instant 
the walker jumps a distance 2* left or right irrespectively 
with probability 1/2*+^, and we can write 

'Pn,t ^ 

(87) 

with 

^ oo 

Tn,n' = r,^ ' ((^Ti-ri'.b' + <5ri-n',-6') ' (^^) 

i=0 

where a = b = 2. Eqs. (j87ll88p are identical to the Weier- 
strass random walk discussed in Refs. [isL [l9l | for ar- 
bitrary 1 < a < b^. There, it was shown that dw = 
ln(a)/ ln(6), which leads to the conclusion that = 1 in 
Eq. pT)) for HN4, as has been predicted (with logarith- 
mic corrections) on the basis of numerical simulations in 
Ref. 



V. CONCLUSIONS 

We have show how the powerful tools of the dynamic 
renormalization group [l^ allow to dissect this intricate 
random walk problem on the planar network HN3 with 
a "hidden" fixed point. Indeed, using a boundary-layer 
analysis, we unravel the irregular singularity of the dom- 
inant fixed point in a five-dimensional parameter space, 
resulting in a set of exact, non-trivial exponents describ- 
ing super-diffusive transport. Adding just one more link 
to each site, wc obtain a non-planar network HN4 which 
possess an even higher degree of symmetry, yet, for which 
we can only develop an equation for the generator and an 
alternative "annealed" treatment which provides results 
that are consistent with simulations. (We believe that 
a proper exploitation of the symmetry in HN4, which 
eludes us here, will ultimately make exact results possi- 
ble.) 

Aside from the singular fixed point, HNS serves further 
as an instructive example for a network in which nodes 
have heterogeneous recurrence properties. The diffusion 
exponent d^ is larger than the fractal dimension d/ = 1 
of the lattice backbone that the walk is embedded in, 
which usually implies recurrence [23,[2l|. Here, the near- 
confined state of the walk favors recurrences to sites in 
higher levels of the hierarchy, although the associated 
first-return exponent is the same for all sites for the time 
distribution of any given return. 

We should also mention that our results for HNS can 
have an alternative interpretation. If we ignore the one- 
dimensional lattice backbone and instead consider the 
network as graph without particular embedding, then 
Eq. ^ for the diameter, or more specifically the aver- 
age growth in neighborhood Sd ^ d^ with jump-distance 
d found in Fig. [51 implies that the fractal dimension 
for that graph is d/ = 2. The RG would discover 
the then-obscured asymmetry between the backbone and 
long-range jumps (even when starting with p = 1/S) 
and lead to the same analysis. Yet, with all distances 
now being (on average) measured as the square-root 
of their separation along the backbone, also the mean- 
square displacement in Eq. (|17|) needs to be reevaluated, 
yielding a diffusion exponent twice its previous value, 
dw = 2(2 — log2 <j)) = 2.61 .... In this interpretation, 
dw > df still applies, but walks are now sub-diffusive in 
this measure. Of course, an exponent that is independent 
of such a metric, like the purely event-base first return 
probability, does not change. In turn, the relation be- 
tween dw and /i fails, consistent with the fact that the 
walk can no longer be considered a Levy flight. 

Finally, our results suggest that many other interest- 
ing transport phenomena, such as voter models, exclu- 
sion processes, or self-organized critical phenomena can 
be fruitfully studied on these networks, which are suf- 
ficiently complex for interesting results but sufficiently 
simple to be tractable. Especially in light of the tremen- 
dous interest in complex dynamics on designed struc- 
tures, we hope that these networks can make a useful 
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contribution 0, S, i, HI . 
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